ApJAccepted 



Normalization of the Matter Power Spectrum via Higher-Order 
Angular Correlations of Luminous Red Galaxies 

Ashley J. Ross 1 , Robert J. Brunner 1 ' 2 , Adam D. Myers 1 
aross2@uiuc . edu 

ABSTRACT 

We present a novel technique to measure a s , by measuring the dependence of 
the second-order bias of a density field on as using two separate techniques. Each 
technique employs area-averaged angular correlation functions (u>n), one relying 
on the shape of a} 2 , the other relying on the amplitude of s 3 (s 3 = oj^/uj^)- We 
confirm the validity of the method by testing it on a mock catalog drawn from 
Millennium Simulation data and finding a ^ easured - a% ue = -0.002 ± 0.062. We 
create a catalog of photometrically selected LRGs from SDSS DR5 and separate 
it into three distinct data sets by photometric redshift, with median redshifts 
of 0.47, 0.53, and 0.61. Measurements of c 2 , and a$ are made for each data 
set, assuming flat geometry and WMAP3 best-fit priors on Q m , h, and T. We 
find, with increasing redshfit, c 2 = 0.09 ± 0.04, 0.09 ± 0.05, and 0.09 ± 0.03 
and a 8 = 0.78 ± 0.08, 0.80 ± 0.09, and 0.80 ± 0.09. We combine these three 
consistent as measurements to produce the result a 8 = 0.79 ±0.05. Allowing the 
parameters Q m , h, and T to vary within their WMAP3 la error, we find that 
the best-fit erg does not change by more than 8% and we are thus confident our 
measurement is accurate to within 10%. We anticipate that future surveys, such 
as Pan-STARRS, DES, and LSST, will be able to employ this method to measure 
erg to great precision, and will serve as an important check, complementary, on 
the values determined via more established methods. 
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Introduction 



The normalization of the matter power spectrum is parameterized as the rms mass 
fluctuation within a top-hat radius of 8/i _1 Mpc and denoted <r 8 . Measuring the amplitudes 
of the matter power spectrum, and thus <7g, is complicated by the fact that most of the 
matter in the Universe is dark and that we must therefore rely on "tracers" of the matter - 
in most cases galaxies. There is no guarantee, of course, that galaxies will cluster in the same 
manner as dark matter. The rela tionship betw een the clustering of galaxies and dark matter 



is known as the "bias" (see, e.g.. iKaiserl Il984l ) . The bias essentially shifts the amplitudes of 



the galaxy power spectrum relative to the matter power spectrum and there is thus a strong 
degeneracy between the bias and a%. Therefore, precisely determining <7g is important, as 
until this is accomplished, the full form of the relationship between the clustering of dark 
matter and the clustering of galaxies will remain ambiguous. 

Measurements of <7g made using cluster counting techniques find mixed results. Using 
the X-ray temp erature and luminosit y functions and fully marg inalizing over the cluster 
scaling relation. IPierpaoli et al.l (120031 ) found <j 8 = 0.77lo!o4> while iHenryl (120041 ) found erg = 
0.66 ± 0.16. Using a self-calibration tech nique and the r ed-seq uence to optically identify 
clusters in the red -sequence cluster survey, bladders et all (120071 ) found a H = 0.67jffjg, while 
Rozo et all J2007h found a 8 = 0.92 ± 0.10 using the SDSS maxBCG fcoester et aliboOTT ) 
catalog. Clearly, there is large variation in the measurements of <7g determined via cluster 
abundances. 

Measurements made using data from the Wilkinson Microwave Anisotropy Probe (WMAP) 
place some of the best constraints on the value of <r 8 , yet still all ow a wide range of values. 
The first year WMAP (WMAP1) results found a 8 = 0.92 ±0.10 JSpereelet, al~1l2003h . while 
the third year (WMAP3) results determined 0.744±g;g| JSpergel et aliboOTT ). The best-fit 
WMAP3 results vary significantly depending on the adopted constraints and priors. In the 
currently accepted cosmological paradigm, inflation is a key ingredient. When WMAP3 is 
constrai ned by a variety of in flationary models, the best-fit a$ are found to be as low as 0.702 
±0.062 (ISpergel et al.l 120071 ) . If one instead combin es the WMAP3 resu lts with the Sloan 



Digital Sky Survey (SDSS) galaxy power spectrum (ITegmark et al.l 120041 ). cig = 0. 772^042- 



Considering the importance of inflation to our understanding of the Universe, the constraints 
placed by WMAP3 are quite loose. Analysis of the WMAP 5 year data (WMAP5; made 
public d uring the revision pro cess of this work) yielded a best-fit 5-year mean value of 0.796 
± 0.036 (IKomatsu et al.l 120081 ) . Despite the precision of this value, if one considers the high 
and low values measured by WMAP1 and WMAP3 and the range in values determined via 
cluster counting techniques, the true value of <t 8 remains unclear. 



In this paper, we present a novel technique to measure <r 8 . Our technique relies on the 
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fact that the bias relationship may be non-linear and if the bias is non-linear, then both 
the amplitude and the shape of the correlation function (the Fourier transform of the power 
spectrum) are affected. As the value of erg affects only the amplitude of the correlation 
function measurement, the non-linear bias is not degenerate with c 8 . Therefore, different 
techniques for measuring the extent of the non-linear bias can constrain not only the non- 
linear bias, but erg as well. 

The bias relationship can be expressed as a Taylor expansion, with the parameters bi and 
b 2 representing the firs t- and second-ord er contributions (and b 2 thus being a measure of the 
non-linearity; see, e.g.. lGaztanagalll992l ). In general, it is convenient to express the second- 
order contribution as c 2 = b 2 jb\. The average overdensity increases as the scale becomes 
smaller, therefore the effect of c 2 increases on smaller scales. If c 2 is positive, the correlation 
function will be increasingly amplified towards smaller scales. Thus, transitioning between 
the non-linear and linear regime (~10 h Mpc, the "weakly non-linear regime") the bias 
relationship changes such that b 2 transitions from having little effect on the measurement 
to having an important effect. Therefore, it is ideal to test c 2 using correlation functions 
at scales on either side of the weakly non-linear regime, as the effects of b\ and b 2 can be 
decoupled by considering both small and large scales. 

Photometrically selected luminous red galaxies (LRGs) from the SDSS DR5 imaging 
data are ideal for measuring correlation functions in the weakly non-linear regime. The 
median redshift of SDSS LRGs is ~0.52, meaning that 10 ft. -1 Mpc is equivalent to about 0.4 
degrees — a scale at which angular galaxy correlation functions in SDSS can be calculated 
quite precisely. Thus, by studying the shape of the 2-point area-averaged correlation function 
{uj 2 ) of LRGs around these scales, one can measure c 2 quite precisely. This measurement is 
dependent on erg and therefore one can determine the relationship between the non-linear 
bias of LRGs and erg. A separate relationship between the non-linear bias and erg can be 



found by utilizing the shape of the hierarchical amplitude s 3 = u 3 / {uj 2 ) 2 (see, e.g.. lRoss et al. 



20071 ; R07 from hereon). These two relationships are nearly orthogonal, thereby allowing 
unambiguous determinations of b\, c 2 , and erg. 



This i s not the first study of the clustering of photometrically selected LRGs. iPadmanabhan et al 



( 120071 ) and lBlake et al.l (120071 ) selected LRGs from SDSS imaging data and studied their clus- 
tering. These efforts measured the two-point correlation function and power-spectrum, re- 
spectively, at different redshifts to place tight constraints on the matter and baryon densities 
of the universe ( but not erg). No r is th is the first study to measure the higher-order cluster- 
ing of LRGs, as iKulkarni et al.l (120071 ) measured the redshift space three-point correlation 
function of spectroscopically selected LRGs. 



Our work, however, has three major distinctions from these previous studies: (1) our 
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data is drawn from the contiguous area of SDSS DR5, significantly improving our ability to 
calculate higher-order statistics at large scales, (2) we measure higher-order statistics using 
photometrically selected LRGs and thus the median redshift is higher than for any previous 
sample used to study the higher-order clustering of LRGs, and (3) we focus on determining 
the value of as- We will outline our methods for determining the two separate relationships 
between as and C2 in §2. In §3, we will describe the creation of a photometric redshift 
ca talog of LRGs d rawn from the SDSS DR5 imaging data, which follows the prescriptions 
of ICollister et al.l (120071 ) . In §4, we will present precise measurements of the second- and 
third-order angular area-averaged correlation functions of LRGs. We will display our mea- 
surements of the first- and second-order bias parameters of LRGs, and we will show that 
these bias parameters have a strong dependence on as- We combine these measurements to 
present a precise measurement of cig. 

Unless otherwise noted, we assume a flat cosmology with parameters equal to the 
WMAP3 alone best-fit values (ft m , h, T) = (0.238, 0.73, 0.135), where T is the shape pa- 
rameter, and essentially parame terizes assuming fixed Vt m and h (see, e.g., Equations 30 
and 31 of lEisenstein fc HiJIlQQsh . 



2. Methods 

2.1. Angular Correlation Functions 

We estimate N-point area-averaged angular correlation functions, ijn(9), using a counts- 
in-cells technique identical to that used in R07. This involves calculating the statistical 
moments of the over-densities contained in equal-area cells. We create the cells using a 
modified version of theSDSSpix pixelization scheme originally developed by Tegmark, Xu, 
and ScrantoK0. The over-density for cell i is defined as 

Si = ^ (1) 
n 

where n is the average number of galaxies in a cell and is the number of gal axies in cell 



i. The remaining details and equations required to determine uj n {6) are found in lRoss et al. 



(120061 ) . The hierarchical amplitudes are defined as 



s * = -J^T I 2 ) 



1 http:/ /lahmu. phyast.pitt.edu/~scranton/SDSSPix/ 
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2.2. Bias and cr 8 



In order to determine the relationship between the bias parameter c 2 and erg, measure- 
ments of uj 2 and S3 are compared to theoretical models constructed via matter power spectra. 
In each case, the z = value of <r 8 is input to the model, allowing c 2 (cr 8 ) to be calculated. 



Model u 2 are produced using power spectra calculated using the ISmith et al.l (120031 ) 



fitting formul ae, as described in R07 §6. By using a modified version of Limber's equation 



(jLimberlll954j ). one can use the redshift distribution of our LRG catalog to invert the P(k) 



and obtain u) 2 (Q): 

u 2 (6) = ^- J J ^ + zf + Q A P(k, z)W 2 2 D [ X (z)9k}dzdk (3) 

where W 2 d — 2'-^^- is the top-hat two-dimensional window function, x( z ) 1S the comoving 
distance to redshift z, P(k, z) is the matter power spectra, k is the spectral index, and J\ 
is th e first-order Bessel func tion of the first kind, and this equation requires Qtotai — 1 (see, 
e.g., Bernardeau et al. 20021 ). This equation is then integrated using the assumed cosmology 



and the desired value of cr 8 . 

Both methods employed to measure c 2 are dependent on the product of bi and <r 8 . In 
order to account for this, we determine the first-order bias using the ui 2 measurement and 
the model u 2 at cr 8 = 0.8 and denote it &i,o.8- This is calculated for scales where linear theory 
is a good approximation (> lOh^ 1 Mpc). Therefore, a valid expression for the first-order 
bias is given by b liLRG = (0.8/cx 8 ) & 1>0 . 8 . 

The second-order bias can be determined by manipulating the overdensities in each cell 
used in the calculation of u) 2 . The overdensity of LRGs can be related to the overdensity of 
dark matter and bias terms via a Taylor expansion (to second order) 

Slrg = bi5 DM + 0.5b 2 5 2 DM (4) 

In order to apply a second order bias term to the LRG correlation measurement, one must 
solve Equation H] for 5dm, which, to second order, is 

8dm = hRc/h - 2b 2 b\ RG jb\ (5) 



Thus, in order to determine the first and second order bias of the LRGs, we use Equation 
Oto apply a b\ and a b 2 to each overdensity used in the measurement of Cu 2 and match these 
altered measurements to the model ui 2 . To find the best fit b 2 for a given <r 8 , one must 
simply calculate uj 2 for a sampling of b 2 values, calculate the \ 2 for each, and minimize x 2 
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via iteration (fully accounting for covariance as noted in §4). This must then be repeated 
for all of the <Tg values one wishes to test (which requires determining the model u 2 for each 
cr 8 ). This process requires the correlation functions to be calculated thousands of times, but 
is highly parallel. In practice, we determine the \ 2 values for selected a%/b 2 values on a grid 
of initial spacing 0.02 in ag and 0.005 in b 2 . This grid is then refined in areas of rapidly 
changing \ 2 ■ The x 2 value for any as/b 2 pair is then found by using a two dimensional spline 
fit to the grid. This method will hereon be referred to as the shape method. 

The other method w e employ is detailed in R07, and takes advantage of the relationship 



f lFrv fc Gaztanagalll993h 



S3 = W (s 3jDA f + 3c 2 ) (6) 

where S3 is the measured amplitude and s 3 ,dm is the theoretical amplitude. The s^dm are 
determined at scales greater than 8 /i -1 Mpc. This is accomplished by calculating ) 3 via the 



integration of linear power spectra and redshift distributions given by iBernardeaul (119951 ). 

Us,dm = 6 (^f ) 2 / (^) 3 [fUl + + ^a] dz X 

{§ (/ kP(k)W$ D [D6k]dk) 2 + J kP{k) {Wi D [D6k}) dk j k 2 D6P{k)W 2D [D6k]W 2D [D6k]dk} 

(7) 

where D is the comoving distance to the median redshift. We calculate uj 2 using Equation [3] 
and linear power spectra and then use s 3 = u)^jQ\. By altering <r 8 and using the relationship 
described in §7.1, Equation 14, of R07, we can determine the la allowed region of c 2 /as- We 
will refer to this approach as the R07 method. 



2.3. Testing via Mock Catalogs 



To test our new method, we took galaxies wit h M r < -23 from th e Blaizot all-sky 
mock catalog cre ated using the metho ds described in iBlaizot et al.l (120051 ) and Millennium 
Simulation data (ISpringel et al.l 120051 ) . In order to select red galaxies we constrained the 
absolute magnitudes of our mock catalog to have B — R > 1.4, as we discovered that the 
color of the simulated galaxies was bimodal about this value. This yielded a sample of 
nearly 300,000 simulated LRGs. These simulated LRGs had a median redshift of 0.2; this is 
significantly smaller than the LRGs to be used in our measurements, but nonetheless quite 
sufficient to test our measurement techniques. 



Using our mock catalog, ujn and sn were calculated using the methods described in §2.11 
In order to determine the bias of the simulated LRGs, we employed the methods described 
in §2, to calculate model u>n and s^- This was done using cr 8 = 0.8 and the assumed 
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cosm ology of the Millenni um Simulation (relevant parameters being fl m , h, T — 0.25, 0.73, 



0.14; ISpringel et al.ll2005f ). This allowed us to find & 1)0 . 8 = 2.04 ± 0.02, fit at scales > 8.2 
h~ 1 Mpc (> 0.8°). With this value in hand, we could then use the shape method to find the 
X 2 values in the b 2 /a 8 parameter space. This allowed us to find the best-fit c 2 value as a 
function of erg. These la bounds, which represent A% 2 = 1 from the minimum at that erg, 
are displayed on Figured] by solid black lines. For <jg = 0.8, we found c 2 = 0.186 ± 0.026. 

Using the R07 method and setting a$ = 0.8, we found c 2 = 0.41 ± 0.09. The large 
disagreement with the value of c 2 determined by the shape and R07 methods is expected for 
a$ = 0.8, as the two methods should agree only for the as used to create the mock catalog 
(0.9). To find where the methods agreed, we calculated \ 2 f° r t ne entire c 2 /<Tg parameter 
space using the R07 method. The resulting la allowed region, which represent Ax 2 = 1, is 
bounded by the dashed black lines in Figure [TJ Combining the \ 2 distributions of the two 
methods produced the la (Ax 2 = 2.3 from the overall minimum) region displayed in red in 
Figure[TJ From these measurements, we determine erg = 0.898 ±0.062 and c 2 = 0.146 ±0.037. 
This measured as is entirely consistent with the input as = 0.9 of the Millennium Simulation. 
This confirms that our method can indeed be used to measure the value of <jg both precisely 
and accurately. 



3. Data 



We take data from the fifth data release (DR5) of the Sloan Digital Sky Survey. To 
cre ate a catalog o f LRGs with photometric redshifts we applied the techniques described 
by ICollister et al.l (120071 ) (abbreviated C07 from here on) applied to objects in the DR5 
PhotoPrimary view. Employing the color and magnitude cuts described by C07 produces 
a sample of just over 1. 7 million objects . As in C07, we found photometric redshifts by 
using the annz software (IFirth et al.ll2003T) with the Two-D egree Field SDSS LRG and QSO 
(2SLAQ) spectroscopic LRG catalog (jCannon et al.l 120061 ) . with stars removed, as training 
data. In order to separate stars, we again used the annz software, and trained it on the 
2SLAQ LRG target catalog. In this case, we included the targeted objects determined to 
be stars and gave them a classification value of 0, while galaxies were given a classification 
value of 1. This same method was employed by C07 to eliminate stars from their catalog. 

Our final catalog comprises only objects with classification values greater than 0.8. 
Based on the training data, cutting on this value should reduce stellar contamination to less 
than 2% while keeping 99.9% of the LRGs. This results in a catalog of 1,662,390 LRGs 
with a median photometric redshift of 0.52. Our redshift distribution is nearly identical to 
the distribution found by C07. These LRGs are then processed through the same imag- 
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ing/reddening/seeing masks as in R07, leaving 1,168,702 objects. We split these LRGs into 
three distinct photometric redshift ranges with similar numbers of objects, 0.4 < z < 0.5 
(444,175 LRGs), 0.5 < z < 0.57 (398,250 LRGs), and 0.57 < z < 0.7 (326,277 LRGs). 
These data sets will be referred to as Z QA1) Z 53 , and Z 061 , with median redshifts of 0.47, 
0.53, and 0.61, respectively. This gives us three distinct data sets that allow us to test the 
consistency of our measurements and that can be combined to increase the precision of our 
final <Tg measurement. 



4. Measurements 



We calculate the area-averaged angular correlation functions (ujn) and hierarchical am- 
plitudes (sjv) for photometrically classified SDSS DR5 LRGs using the methods described 
in §2.11 For every measurement, errors an d covariance matrices are calculated using a jack- 



knife method 
(e.g., 



Scranton et al.l 120021 ) , with inv erse-variance weig hting for both errors 



Myers et al.ll2005l 120061 ) and covariance (e.g., iMyers et al.l 120071 ) . identical to the one 
described in §3.4 of R07. This allows us to minimize x 2 i fully accounting for covariance via 



X 



5>(0i) - M0*)]ctf[0(*i) " 



1,3 



where C is the covariance matrix, and i and j refer to the i and j jackknife subsample. 

Focusing first on Z 053 , we fit u 2 for measurements made between 0.4° and 1.6° (10.2 
/i _1 Mpc and 40.2 ftr^Mpc). We determine & 1(0 . 8 = 1.63 ± 0.02. The Q 2 is well fit by a 
single bias parameter in this range, as the x 2 = 1-7, P{< X 2 ) = 0.89. Thus, we measure 
bi,LRG = 0.8/o"8 x (1.63 ± 0.02). The u)i measurement is presented in Figure [21 divided 
by 1.63 2 (accounting for fr^o.s) along with the model a> 2 at <r 8 = 0.8. At scales less than 
0.3°, the measurement grows larger than the model, indicative of positive second-order bias. 
The measurement also grows larger at scales greater than ~2° (50.2 /i -1 Mpc), but at these 
scales the errors begin to grow larger and systematic s due to reddening and projection effects 
also increase. Based on the results of ISimonl (120071 ). our theoretical curve, which employs a 
modified version of Limber's Equation, should not be accurate to better than 10% at scales 
greater than ~2° for any of the redshift ranges we use. Our measurement at 2° diffe rs from 
the m odel by 8.5%, thus the disagreement is no greater than would be predicted by [Simon 
( 120071 ). We thus fit no measurements to scales greater than 1.6°. 



Altering the C02 measurement using the shape method, and fitting the measurements at 
scales between 0.1° and 0.7° (equivalent to ~ 2.5 to 17.6 h^ 1 Mpc, there are 8 measurements 
in this range and thus 7 degrees of freedom), we measure 62 = 0.15 ± 0.05 for a$ = 0.8. The 
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fit is acceptable, as \ 2 — 0.60, P{< y 2 ) = 0.999. Attempting to fit the data with a single 
bias parameter model, we find x 2 = 12.18, P(< x 2 ) = 0.09. For other redshift ranges, we 
find similar results. Again for a 8 = 0.8, we measure &2,z .4 7 = 0.150 ± 0.040 (x 2 = 0.49, 
X 2 ) = 0.999). A single bias parameter model is rejected at 87%. Finally, we find 
^2,z . 6 i = 0.165 ± 0.025 (x 2 — 1-91, P(< x 2 ) — 0.96). A single bias parameter is rejected at 
> 99% for this redshift range. Based on the marginal rejections of a single-parameter model, 
a two-parameter model is needed to fit the measurements for each redshift range. For the 
two lower redshift ranges, the minimum x 2 values are quite small, implying that perhaps our 
error bars are overestimated for these redshift ranges, further implying our quoted errors on 
6 2 are overestimated. 

In the bottom right panel of Figure [31 the u>2 measurement corrected for b\ = 1.63 and 
62 = 0.15 is displayed along with a theoretical for cr 8 = 0.8. The model curve clearly 
fits the data. The other panels display the measured S3 (black triangles) for Z 0A7 , Z , 53 , 
and Z .6i (left to right, top to bottom) corrected for the best-fit b\ and C2 in accordance 
with Equation [6j Each panel also includes a solid line displaying the model S3. The Z0.53 
measurement appears to be extremely consistent with the model, while the curves defined 
by the other two measurements do not share the same shape as the model. Despite this fact, 
the size of the error bars allows the Z0.47 and Z0.61 measurements to appear consistent with 
the model. 

We use the R07 method to find c 2 for each data set. For data sets Z 0A7 and Z 053 we 
fit between 0.4° and 1.6° (equivalent to 9.0 h^Mpc and 35.9 h^Mpc for Z 0A7 and 10.0 
/i _1 Mpc to 40.2 /i _1 Mpc for Z0.53; there are 7 measurements in this range and thus 6 degrees 
of freedom). For Z0.61, we fit between 0.3° and 1.6° (8.5 /i _1 Mpc to 45.6 /i _1 Mpc; again 7 
measurements and 6 degrees of freedom). We find that for cr 8 = 0.8, c 2) z 047 — 0.08 ± 0.12, 
c 2 Zo 53 = 0.07 ± 0.13 and c 2i z . 6 i = 0.14 ± 0.16. For Z . 53 , x 2 = 0.033, meaning P{< X 2 ) = 
1.0 — 1.0 x 10~ 6 . This is a remarkably small x 2 value, which one might expect (to a degree) 
based on how well the measured values appear to match the model, despite the size of the 
error-bars. This suggests that the error-bars are over-estimated for the S3 measurements in 
this redshift range, and thus the error on our C2 measurements may be overestimated as well. 
For both the Z 0A7 and Z a61 ranges, x 2 — 1-5 with P{< x 2 ) — 0.96. These values are quite 
reasonable, implying that if our errors are being overestimated, it is happening only for the 
Z0.53 range. 

Using the shape and i?07methods (see §2.21) the x 2 values over the entire cr 8 /c2 parameter 
space are determined for each data sample. The la allowed regions of c 2 /cr 8 determined via 
the R07 (solid black lines, A% 2 = 1 from fixed <r 8 ) and the shape (dashed black lines, A% 2 = 1 
from fixed <r 8 ) methods are plotted in Figure HJfor Zq.53- Fortunately, the two methods bound 
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significantly different regions of parameter space, allowing a precise determination of c 2 and 
erg. Combining the two measurements produces the la (red triangles, A% 2 = 2.3 from overall 
minimum) allowed regions for c 2 /as, and are also displayed in Figure HI To la precision, we 
thus find c 2j z 53 = 0.092 ± 0.052 and a 8 = 0.796 ± 0.086. Repeating the process for Z 0A7 , we 
find c 2 ,z 0A7 = 0.088±0.041 and a 8 = 0.776±0.080 and for Z .ei, we find c 2 ,z . 61 = 0.092±0.033 
and a 8 = 0.798 ± 0.094. 

The best fit <7g for our three data sets are consistent to 0.275cr. Combining the three 
measurements, we measure erg = 0.789 ± 0.050. Adopting this value in order to determine 
the first order bias, we find &i,z 047 = 1-47 ± 0.09, &i,z 053 = 1-65 ± 0.09, and &i,z 61 = 
1.80 ± 0.10. These values make sense given that the median luminosity of the galaxies 
increases with redshift since our sample is not volume limited. If we multiply each of the 
best-fit c 2 measurements by 0.789 divided by the best-fit erg for each respective data set 
(approximately correct for small changes in erg, based on our shape method measurements), 
we find that c 2 , ZoA7 = 0.09 ± 0.04, c 2 ,z . 5S = 0.09 ± 0.05, and c 2 ,z . 61 = 0.09 ± 0.03 for 
ag = 0.789. The fact that there is no significant change in c 2 is moderately surprising and 
implies differences in the halo occupation distribution (HOD). (I f the HOD was not cha nging 



as a function of halo mass, c 2 would increase with see, e.g., iNishimichi et all 120071 ). We 
will discuss the HOD more in ^5.21 



Of interest is the fact that the S3 measurement for both Z0.47 and Z0.61 have a local 
minimum (displayed in the two left panels of Figure [31 but there is no such minimum in the 
S3 of Z 053 . The minimum is at ~0.6°, equivalent to 11.4 /i _1 Mpc for the lower redshift range, 
while it is at ~0.3°, equivalent to 8.5 /i _1 Mpc, in the high redshift range. In R07, it was 
found that early-type galaxies also displayed a minimum in their S3 measurement at ~10 
h^Mpc. As in R07, the errors dominate the LRG measurement (though to a lesser extent), 
but it appears unlikely that this is a coincidence. The feature is seen at approximately 
the same physical scale, but at a different angular scale due to the differences in redshift. 
This rules out any possibility of observational systematics such as seeing or reddening. It is 
unclear whether the feature may be due to projection effects or complicated halo dynamics 
that do not affect the middle redshift range. 



5. Discussion 

We have presented a technique for measuring as using the 2 and 3-point angular area 
averaged correlation functions and applied it to photometrically classified LRGs from the 
SDSS DR5, split into three distinct redshift ranges. Using a method that depends on the 
shape of uj 2 and the technique described in R07, two separate relationships were determined 
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between a 8 and C2,lrg- These relationships split the degeneracy between bias and a 8 , allow- 
ing an unambiguous determination of the 1st and 2nd order bias and a 8 . The measured a 8 
in three redshift ranges are consistent and combine for a best-fit a 8 = 0.789 ± 0.048. Our 
determination of a 8 is quite precise. It is thus important to investigate the assumptions, 
implicit and explicit, made when determining our measurements, to compare our measure- 
ments to the relevant dark- matter-halo/bias theory, and to investigate how consistent our 
measurements are with previous results. 



5.1. Assumptions 

The main assumption that goes into our measurements of c 2 is that the bias can be 
expressed solely as a function of the overdensity, i.e. it is not a function of both the overden- 
sity and the smoothing scale. If the bias was a strong function of scale, it would invalidate 
any measurement made using the shape method, as a changing bi would change the shape 
of £D 2 . Similarly, the R07 technique assumes a constant bi and c 2 over the range of scales 
that are fit. Further, the R07 measurement is fit at a different range of scales than the shape 
measurement. Bias that is a strong function of scale between 2.5 /i^Mpc and 40 h^Mpc 
would completely invalidate any comparison between measurements using the shape and the 
R07 method, and thus invalidate our a 8 measurements. 

Our assumption that the bias can be expressed solely as a function of the overdensity is 
validated by the goodness of the fit to our bias-corrected measurements. These measurements 
demonstrated that, between 0.1° and 1.6°, our lj 2 measurement is quite consistent with a two 
bias parameter model, as all of the model fits using the shape method were accepted to better 
than 96%. Further, in each redshift range, the probability that a single bias parameter fits 
the data is less than 9%. The simplest model that fits the data between 0.1° and 1.6° is thus 
that the bias is independent of scale and can be described by two parameters. Therefore, we 
believe our measurement techniques to be valid and that our comparison of those techniques 
is valid as well. 

The error on our a 8 measurement is quite low in part because we held the values of Q m , 
h, and V fixed to their WMAP3 best-fit values. Allowing these values to change does alter our 
best fit values of c 2 and <t 8 . In order to determine the degree to which the uncertainty in these 
parameters should affect the uncertainty of our results, we repeated our a 8 measurements 
using the Z0.53 data set (0.5 < z < 0.57) and produced model u>2 and S3 for each parameter at 
the la limits determined by the WMAP3 alone best-fits (while holding the other parameters 
at their best-fit values). 
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As long as the geometry of the Universe is kept flat, changes in the matter density have 
little effect on either the amplitude or the shape of uj 2 . Thus, we expected that the value of 
Q m would have little effect on our measurement of a 8 . This was indeed the case, as we found 
<7 8 = 0.796±0.086 and 0.800±0.088 for tt m equaling 0.251 and 0.2134, respectively, with n total 
fixed at 1 and (T, h) = (0.135, 0.73). Conversely, we expected our measurement of ag to be 
fairly dependent on the value of h, as this value significantly affects the distance to the LRGs. 
We found a 8 = 0.846±0.088 and 0.748±0.085 for h = 0.7 and 0.76, respectively, with {Vt totah 
T, Q m ) = (1, 0.135, 0.238). The percentage change in the measured a s was approximately 
the same as the percentage change in h, suggesting a close relationship between the two 
(again, as expected). 

We also expected the value of T to have a significant effect on the measurement, as 
changing T alters the shape of u> 2 . For T = 0.149 and 0.12, we found a 8 = 0.842 ± 0.088 and 
0.736±0.082, with (tt totah h, ft m ) = (1, 0.73, 0.238). While the uncertainty in T produced the 
largest range in a 8 values, the percentage change in a 8 was actually ~25% smaller than the 
percentage change in T. We thus determine that although the uncertainty in T introduces the 
most uncertainty into our measurement of a 8 , our measurement technique is most sensitive 
to the value of h. Despite the changes in the value of a$ that we measure, our measurement 
of a 8 has not changed by more than 8%. This suggests that our quoted uncertainty of 0.05 
would increase by less than a factor of two when uncertainties in h and T were taken into 
account. We are thus confident our measurements are accurate to within 10%. 

We also assumed no error in our redshift distribution when making our measurements. 
Precise knowledge of the redshift distribution is necessary for Equations 3 and 7. In order to 
explicitly test our measurements' dependence on the redshift distribution, we created two new 
distributions for the Z 52 redshift range (see R07 for the details of how these distributions 
are constructed). For one, we systematically increased the error of each photometric redshift 
by 10% (effectively broadening the redshift distribution) and for the other we decreased the 
error of each photometric redshift by 10% (effectively narrowing the distribution). We then 
re-calculated the best-fit a 8 for each distribution, finding that it increased to 0.84 ± 0.09 
for the distribution with greater redshift errors and that it decreased to 0.76 ± 0.09 for the 
distribution with smaller redshift errors. These differences are significant, but are smaller 
than our combined la error. We thus do not believe this issue significantly adds to our 
quoted uncertainty. 
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5.2. Testing via Halo Models 

It is important to determine how our measured bias values compare to t heoretical values. 



We ca lculate b\ and C2 using halo models and the methods prescribed by iNishimichi et al. 



(120071 ). They show that one can combine the iVth-order bias coefficient of halos as a func- 
tion of mass and redshift, B^(m,z), the number of halos of a certain mass and redshift, 
nhaio(fn, z), and the mean number of galaxies occupying a halo of a certain mass, (n g \M), 
to find the bias of a population of galaxies. This can be expressed by 



>N 



J dMn halo (M,z)B N (M,z) (n g \M) 
j dMn halo (M, z) 



(9) 



We determine b oth n ha i (m,z) and B N (m,z), by using an ellipsoi dal collapse mod e l (e.g ., 



Sheth et al.ll200ll ) and following the methods described in detail in INishimichi et al.l (120071 ). 
For (n g \M) we first model the number of central LRGs per halo as having a "soft" transition 



between n n = and n„ = 1 such that 



(ncentral\M) =0.5 



1 + erf 



\og w {M/M cut ) 



&cut 



(10) 



as m 



Zheng et al. 


(2005 


) and 


Blake et al. 


(2007a) 



are not using a volume-limited sample, and t hus at smaller redsh ifts the mass limit is likely 
to be smaller than at higher redshifts. As in iBlake et al.l (j2007al ). we model the number of 
satellite galaxies using a simple power law. Thus 



<",,<|M> = ( £ 



(11) 



The bias model has four free parameters, M cut , M , a cut , and a. It is beyond the scope 
of this paper to fit for these parameters. It is instead our intention to determine if reasonable 
values for these parameters can reproduce the bias values that we measure. If this is possible, 
it suggests that the our bias measurements are themselves reasonable. Using cr 8 = 0.793, 
z = 0.532 and halo parameters login(M cut /M ) = 13.3, a cut = 0.6, logio(M o /M ) = 14.5, 
and a = 2.0, we find hi = 1.61 and C2 = 0.09 - - results that are consistent with our 
best-fit values. The values of log i n (M n /M ^) and a cut were chosen to be equal to the best-fit 
parameters found bv Blake et al.l (j2007ah for LRGs with 5 < z < 0.55. The logi (M CMt /M Q ) 
and a parameters are slightly lower than the IBlake et al.l (]2007al ) parameters likely because 
our LRGs have a lower luminosity an d thus a lower bia s, mi nimum mass , and a (a has 
been seen to increase with b\ in both IBlake et al.l l2007al and IZehavi et al.l 120051 ) . This is 
not to suggest t hat we favor a steep power-law to a shallow one. If we change a = 1.4 
(as measured by iKulkarni et al.l 120071 ) . we calculate b\ = 1.62 and C2 = 0.09 if we also 
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reduce logio(Mo/M ) to 14.41. This implies that, from a theoretical standpoint and in the 
context of relevant measurements of LRG halo properties, our measurements of hi and c 2 
are reasonable. 

The bias model can also be used to test our implicit assumption that (c 2 ) = (b 2 )/ (bi)- 
We calculate bi, b 2 , and c 2 at redshifts between 0.4 and 0.7 (essentially the redshift range of 
our entire sample) in intervals of 0.02. These values are used to calculate theoretical values 
for (b 2 ) / (b\) and (b 2 /bi). We find that the difference between the two is less that 0.1%, 
meaning that the systematic error introduced by assuming (c 2 ) = (b-z) / (b\) is insignificant 
for our measurements. 



5.3. Comparison With Other Measurements 



Our measurement of <7g = 0.789 ±0.050 is consistent with most previous measurements. 
There are, however, notable exceptions. Our result is inconsistent to la with a photomet- 
ric optical cluster counting technique e mploying SPSS dat a that found a% = 0.92 ± 0.10 
( IRozo et al.ll2007l ). Another recent result (IHarker et al.ll2007l ) found a similarly high value of 
as = 0.97 ± 0.06 using N-body simulations in combination with semi-analytic galaxy forma- 
tion models and the projected two-point correlation function of SDSS galaxies. Other results 
determine <7g to be too small to be consistent with our measurement (to la). Notable exam- 



ples are the W MAP3 constrained by inflationa ry models result of 0.702 ±0.062 (ISpergel et al. 
2007h and the IVoevodkin fc Vikhlinrnl J2004T ) results that found a 8 = 0.72 ± 0.04 by using 
the cluster baryon mass function. 

However, many notable results are consistent with our measurement. The best fit 



WMAP3 data alone found a 8 = 0.744; 



-0.05 
-0.06' 



which is just barely consistent with our mea- 



surement to la. The WMAP3±SDSS and WMAP3±LRG best-fit values of a 8 = 0.772±g;g^| 
and 0.781 

-ao45 are consistent with our measurement to less than la. Large disagreement 
with these results would be surprising, given that we used WMAP3 best-fit priors for the 
re levant input cosmologi cal parameters. The WMAP5 results are quite similar to our results, 



as 



Komatsu et al.l (120081 ) find a$ 



0.796 ±0.036, for the WMAP alone five-year mean value. 
More significantly, our results are also consistent to la with a$ = 0.671^3 as derived from 
optical cluster finding techniques ( Gladders et al.l 120071 ) and a 8 = 0.66 ± 0.16 from x-ray 
cluster measurements ( iHenryl 12004] ) . Due to the fact that the inconsistent measurements 
appear to be as likely to be lower than our measurement as they are higher, we feel that 
these results are hinting at a potential convergence to a a$ value that is close to 0.8. 
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6. Conclusions 

We present and test a new method for determining the value of a§. The method and 
the results of our testing can be summarized as follows: 

• The technique for measuring a 8 utilizes two measures of the second-order bias of a 
density field. The two methods have different dependencies on the value of <t 8 , and can thus 
be combined to determine a best-fit a 8 and second-order bias. One measure of the second- 
order bias (the -R07 method) has been used many times before (e.g., R07) and depends 
on the amplitude of S3. The other method (the shape method) has (to our knowledge) 
never been used before. It relies on correcting the overdensities for given first- and second- 
order bias parameters and determining the bias parameters that allow the shapes of the O2 
measurement and model to become consistent. 

• The method was tested u sing a mock catalog of LRGs drawn from the Blaizot all- 
sky catalog iBlaizot et aL 2005 ) that was constructed using Millennium Simulation data 



( jSpringel et al.l 



2005h and its input cosmological parameters for O m , h, and T. The Millen 



nium Simulation assumed a 8 = 0.9 and we measured as = 0.898 ± 0.062. This measurement 
proved that our method is both accurate and precise. 

• We photometrically selected LRGs from SDSS DR5, deter mined photometric red shifts 



for each LRG, and removed stars following the prescriptions of ICollister et al.l (120071 ). We 
divided this LRG catalog into three samples by redshift, with the separate ranges being 
0.4 < z < 0.5 (Z0.47), 0.5 < z < 0.57 (Z . 53 ), and 0.57 < z < 0.7 (Z . 6 i). 

• We measured a 8 in each sample and found a 8 = 0.776 ± 0.080, a 8 = 0.796 ± 0.086, 
and a 8 = 0.798 ± 0.094, respectively. Combining these consistent results we determined 
erg = 0.789±0.050. For each measurement, we assumed the relevant cosmological parameters 
were equal to their WMAP3 best-fit values. 

• Allowing the relevant cosmological parameters to vary within their WMAP3 la error, 
we found our measurement of a 8 changed by less than 8%. Thus, even allowing for these 
uncertainties, our method produces a precise measurement that we are confident is accurate 
to within 10%. 

• We measured b hZoA7 = 1.47 ±0.09, b 1>Z0m53 = 1.65 ±0.09, and b hZo61 = 1.80 ±0.10 and 
c 2jZoA7 = 0.09 ± 0.04, c 2 ,z .53 = 0.09 ± 0.05, and c 2 ,z . 61 = 0.09 ± 0.03. Using a halo model, 
we determined that the bias values for Z0.53 were consistent with reasonable, and previously 
measured, HOD parameters. 

The techniques described herein can easily be repeated and tested using other cosmic 
samples. Future surveys, such as the Panoramic Survey Telescope & Rapid Response System 
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(Pan-STARRS), the Dark Energy Survey (DES), and the Large Synoptic Survey Telescope 
(LSST), will provide opportunities to measure both the 2- and 3-point area averaged corre- 
lation functions to extreme precision. Repeating the techniques we have presented here will 
allow researchers to determine a 8 accurately, precisely, and independently of future cluster 
counting and CMB techniques, thus providing an important check on those results. 

In addition to our determinations of the bias and of ag, we have found that a feature 
in the hierarchical amplitude of S3 at ~10 /i _1 Mpc exists in two of the three data sets. 
Given that a feature existed at approximately the same physical location for early-type 
galaxies at significantly smaller redshifts (R07), the feature appears physical in nature. The 
fact that it is absent in one of our data sets hints that the possible cause may be due to 
projection effects, or that perhaps the feature is indicative of complicated halo occupation 
statistics. The feature demands further study, both observationally and theoretically, and 
we are currently focusing our efforts to explain this phenomon. 
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Fig. 1. — The la allowed regions for Cija% determined by the shape method (solid black 
lines) and the R07 method (dashe d black lines) for simulated LRGs from the Blaizot all-sky 
mock catalog ( Blaizot et al. 2005 ) created using Millennium simulation data (jSpringel et al. 



20051 ). These measurements are combined to produce the la (red) allowed region. The la 
best-fit values are C2 = 0.15 ± 0.04 and a$ = 0.898 ± 0.062, consistent with the Millennium 
input value of a$ = 0.9 
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Fig. 2. — The area-averaged angular two-point correlation function for photometrically- 
classified LRGs from SDSS DR5 with photometric redshifts 0.5 < z < 0.57 {Z 53 ), divided 
by 1.63 2 (black triangles, fei = 1.6 3 for erg = 0.8). The solid line represents the model u 2 
calculated using ISmith et al.l (120031 ) power spectra and Limber's equation. At scales smaller 
than 0.25°, the measurement is larger than the model, suggesting positive second order bias. 
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Fig. 3. — Left to right, top to bottom, the first three panels display the measured S3 (black 
triangles) for photometrically classified LRGs from SDSS DR5 with photometric redshifts 
0.4 < z < 0.5 (Z0.47), 0.5 < z < 0.57 (Z . 53 ), and 0.57 < z < 0.7 (Z om ). In order to 
take the bias into account, S3 is divided by the best-fit b\ and then has three times the 
best-fit c 2 subtracted from it for each respective redshift range (for a 8 = 0.8). The solid line 
in each of these three panels is the model s 3 for er 8 = 0.8. The Z 53 data has a different 
shape than the other data sets at large scales, but due to the size of the error bars, all data 
sets are consistent with the model s 3 to at least 96%. The lower-right panel displays the 
bias-corrected measurement of cu 2 (black triangles) for Z Q 53 using the best-fit b\ = 1.63, and 
b 2 = 0.15 for <7g = 0.8. The solid black line is the model uj 2 for a 8 = 0.8. Visually, the fit 
appears to be as good as the x 2 — 0.6, P(< x 2 ) = 0.998 suggests. 
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Fig. 4. — The la allowed regions for C2ja% determined by the R07 method (dashed black 
lines) and the shape method (solid black lines) for photometrically selected LRGs from SDSS 
DR5 with 0.5 < z < 0.57 (Z >53 ). These measurements are combined to produce the la (red) 
allowed region. The la best-fit values are C2 = 0.092 ± 0.052 and as = 0.796 ± 0.086. 



